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Recent developments in finite-temperature QCD with dynamical quarks are reviewed focusing on the topics of 
critical temperature, the equation of state, and critical behaviors around the chiral phase transition. 



1. Introduction 

In the last years, much effort has been paid 
to experimentally detect the quark-gluon plasma 
(QGP) in heavy-ion collisions. In these experi- 
ments, precise theoretical inputs about the nature 
of the transition and the equation of state (EOS) 
of QGP are indispensable to unambiguously iden- 
tify the signal of QGP. Here, lattice QCD pro- 
vides us with the most powerful method to com- 
pute thermodynamic properties of QGP from the 
first principles of QCD. 

In the quenched approximation of QCD, the 
study of finite temperature QCD is already well 
matured: Computations of Tc and EOS have 
been made on lattices with several temporal lat- 
tice sizes TVt, using both the standard plaque- 
tte action and various improved actions. Con- 
tinuum extrapolation to the limit of large Nt has 
been performed with different actions to confirm 
good agreement within an accuracy of a few per- 
cent For the pressure, precise comparisons 
of the results from different calculation methods, 
i.e. the integral method[Q and the derivative 
method ||], have also been madej^-^. The ex- 
pected absence of the pressure gap at the first or- 
der deconfinement transition point has also been 
demonstrated through a new non-perturbative 
calculation of anisotropy coefficients 0. 

The quenched study is now entering a new 
stage: Thermal effects on the pole mass have 
been studied by the QCDTARO collaboration us- 
ing an anisotropic lattice ||]. Several groups are 
also trying to compute hadron spectral functions 
at finite-temperature applying the maximum en- 
tropy method||l^-|l^. 

On the other hand, the study of full QCD at 



finite temperatures is still at the stage of de- 
velopment. Until recently, EOS with dynamical 
quarks has been computed only for the staggered 
fermion. The MILC collaboration calculated EOS 
in 2 fiavor QCD with the standard staggered 
fermion on iVt = 4 and 6 lattices . The Biele- 
feld group computed EOS for Nf = 4 using an 
improved staggered fermion at Nt = 4[^. This 
year, the Bielefeld group extended the compu- 
tation to Nf — 2 and 3 and studied the Nf- 
dependencepq]. 

A new development of this year is the first cal- 
culation of EOS for a Wilson- type fermion [pTf. 
Using a renormalization group (RG) improved 
gauge action coupled with 2 flavors of clover- 
improved Wilson quarks the CP-PACS col- 
laboration performed a systematic study of the 
pressure and the energy density of QGP. 

In this report, we concentrate on the develop- 
ment in finite temperature QCD with dynamical 
quarks, focusing on the topics of EOS, and at- 
tempt a comparison of the results from different 
lattice fermions. In section we discuss about 
the phase structure, critical properties around the 
chiral phase transition of 2 fiavor QCD, and the 
transition temperature, we then discuss new cal- 
culations of EOS in section |[ A brief summary 
is given in section ^ 

2. Phase structure 

2.1. Phase structure with Wilson fermions 

It has been difficult to study EOS with Wilson- 
type quarks due to several difficulties: One is the 
complicated phase structure due to the explicit 
violation of chiral symmetry and due to the ap- 
pearance of the parity-fiavor broken phase (Aoki 
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Figure 1. Phase diagram for the RG-improved 
gauge action and the clover-improved Wilson 
quark action on an A^t = 4 lattice p^. 
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Figure 2. 0(4) scaling relation measured on a 
16 X 4 lattice with the RG gauge and clover quark 
actions [E4|. 



phase) |2^j2^ . Therefore, a systematic study sur- 
veying a wide range of the parameter space is 
required. The second reason is that, when the 
standard plaquette gauge action and the standard 
Wilson quark action are used, lattice artifacts are 
large on coarse lattices used in most finite tem- 
perature simulations. For example, unexpected 
strong phase transition is observed at intermedi- 
ate quark masses in 2 flavor QCD|2|,|2|. There- 
fore, we have to use improved actions to suppress 
these lattice artifacts. 

The CP-PACS collaboration studied 2 flavor 
QCD with an RG-improved gauge action com- 
bined with a clover-improved Wilson quark ac- 
tion||2^. The first step is to clarify the phase 
structure for this combination of actions. 

Figure |l] is the result for the phase diagram of 
QCD obtained on an A't = 4 lattice. The solid 
line Kc{T = 0) is the location of the chilal limit 
at zero temperature. The pion mass decreases as 
K is increased from small if, and vanishes on the 
line Kc{T = 0). When Nt is fixed as in most 
simulations, the temperature T = (Nta)~^ be- 
comes higher as f3 is increased. The dashed line 
Kt is the pseudo-critical line for the finite temper- 
ature phase transition for Nt = A. The region to 
the right of Kt (larger /3) is the high temperature 
QGP phase, and that to the left of Kt (smaller /3) 
is the low temperature hadron phase. The cross- 



ing point of the Kc{T — 0) and the Kt lines is the 
chiral phase transition point [2^. 

In the region above the Kc{T = 0) line, a 
parity-flavor symmetry of the Wilson-type quark 
action is broken spontaneously!^,^. At zero 
temperature, the boundary of the parity-flavor 
broken phase (the Kc{T = 0) line) is known 
to form a sharp cusp touching the free massless 
fermion point K = 1/8 at f3 = oo. 

At finite temperatures, the parity-flavor bro- 
ken phase retracts from the large /3 region. The 
boundary of the parity-flavor broken phase at 
T > (the KciT > 0) line) is shown in Fig. | 
for Nt = 4. As shown in this flgure, the chiral 
transition point locates close to the cusp of the 
parity-flavor broken phase, the difference being 
at most of 0{a). This is consistent with the pic- 
ture that the massless pion, i.e. the Goldstone bo- 
son associated with spontaneous chiral symmetry 
breaking, appears only in the cold phase. 

2.2. Scaling study of 2 flavor QCD 

An important step towards clarification of the 
QCD transition in the real world is the study of 
the finite temperature chiral phase transition in 
Nf = 2 QCD. When the chiral transition in Nf = 
2 QCD is second order, the transition is expected 
to be universal to that in a 3-dimensional 0(4) 
spin model: With the identifications M ~ {^'^) 
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for the magnetization, h niqa for the external 
magnetic field and t ~ fi—(3ct for the reduced tem- 
perature, where j3ct is the chiral transition point, 
we expect the same scaling behavior as the 0(4) 
spin model. Therefore, confirmation of the 0(4) 
scaling in QCD provides us with a test of the or- 
der of the chiral transition. 

The 0(4) scaling was first tested with stag- 
gered fermion by Karsch and Laermannp^]. The 
study was extended to a wider range of the quark 
mass and lattice sizes by the JLQCD collabora- 
tion[Q and by the Bielefeld group p7[|. However, 
an agreement of the critical exponents between 
the 0(4) spin model and QCD with 2 flavors of 
staggered fermions was not obtained for Nt = 4. 

For the case of Wilson fermion, Iwasaki et 
a/.p8[| investigated the scaling relation; 



(1) 



using the RG gauge action combined with the 
standard Wilson quark action. They identified 
the subtracted chiral condensate defined by an 
axial Ward-Takahashi identity p9|] 



sub 



2mqaZ^(7r(a;)7r(0)), 



(2) 



as the magnetization of the spin model. Here, 
the quark mass ruq is defined by an axial vec- 
tor Ward-Takahashi identity j2^j3^, and the tree- 
level renormalization coefficient Z — {2KY was 
adopted. They found that the scaling relation (|^) 
is well satisfied with the 0(4) critical exponents 
and the 0(4) scaling function. 

Figure ^ is the new result for the case of 
the clover-improved Wilson quark action coupled 
with the RG gauge action, obtained on a 16^ x 4 
lattice[p4| . The vertical axis is M/h^/^ and the 
horizontal axis is tjh}!^^ , where /3 and 8 are 
the critical exponent obtained in the 0(4) spin 
model The dashed line is the 0(4) scal- 

ing function!^. They fitted the data to the 
scaling function adjusting /3ct and the scales of 
two axes. As seen from Fig. ||, QCD data is 
well described by the 0(4) scaling ansatz also for 
this combination of actions. The best fit gives 
/3ct = 1.469(73) for the chiral transition point, 
with I^DF = 0.82. This result suggests that 
the chiral transition is of second order for Nf = 2. 
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Figure 3. Tpc/mv as a function of mps/mv 
with various actions. 



This year, the MILC collaboration published 
their final results of a study of scaling proper- 
tics in QCD with standard staggered quarks |p3| . 
They found that the QCD data are inconsistent 
with the 0(4) scaling, in accord with the previ- 
ous results for the critical exponents [^ , ^ . At 
the same time, a comparison of the results for 
Nt = 4-12 suggests that the discrepancy with 
the 0(4) scaling function may be removed when 
Nt is sufficiently large. 

A different approach to clarify the origin of the 
unexpected results was proposed by Kogut, Lagae 
and Sinclair |34j. They carried out a simulation 
at precisely zero quark mass by adding the addi- 
tional, irrelevant four-fermi interaction to QCD. 
From their results obtained on A^t = 4 and 6 lat- 
tices, they also suggested importance of simula- 
tions at larger Nt. 

For staggered quarks, a comparison with an 
0(2) spin model is also important because the 
symmetry of the staggered quark action is 0(2) 
at finite lattice spacings. Engels et al. calculated 
the scaling function for the 0(2) spin mo del 
They found that the 0(2) scaling function is quite 
similar to the 0(4) scaling function. This implies 
that the 0(2) scaling also cannot explain the be- 
havior of the staggered quarks at small Nt- 

2.3. Transition temperature 

In Fig. 1^, we summarize the results of the 
pseudo-critical transition temperature Tpc in 
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Table 1 

Chiral transition temperature in physical units. The scale is fixed by ^/a 
{Nf ^ 0) QCD, and by my = 770MeV for Nf = 2 and 3 QCD. 
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Figure 4. Pseudo-critical temperature Tpc in a 
unit of ^/a for 2, 3, 2+1 flavor QCD using the 
p4-improved staggered actionp8|. 



Nf — 2 QCD with various actions; plaque- 
tte gauge combined with the standard staggered 
and Wilson quark p6| , ^ , Symanzik gauge with 
p4-improved staggered quark [ p^j38| , plaquette 
gauge with clover quark Symanik gauge with 
clover quark ||3^, RG gauge with clover quark [p4[, 
and plaquette and RG gauge with domain-wall 
quark|^^. The scale is fixed by the vector meson 
mass my at the simulation point. Filled and open 
symbols are respectively the results on iVt = 4 
and 6 lattices. Diamonds and stars are new; the 
diamonds from the p4-improved staggered quark 
by the Bielefeld group 1 38 1, and the stars from the 
RG gauge with clover quark by the CP-PACS col- 
laboration[Q. We note that the new results are 
slightly larger than the previous results from the 
standard staggered quarks presented by circles. 
The Bielefeld group and the CP-PACS collabo- 



ration also attempted chiral extrapolations. The 
results agree with each other within errors in the 
chiral limit, giving Tc ~ 170 MeV for the chiral 
transition temperature. 

The Bielefeld group also investigated the Nf- 
dependence of Tpc on a 16'^ x 4 lattice[^. Fig- 
ure ^ shows their results for 2 (square) and 3 
(triangle) flavors of degenerate staggered quarks, 
normalized by the string tension. From this fig- 
ure, the dependence on Nf seems to be small for 
^/ = 2-3. 

Table |l| is a summary of the transition temper- 
ature in the chiral limit. The quenched results 
are extrapolated to the continuum limit. How- 
ever, those for full QCD, - 150-180 MeV, are 
the results for Nt = 4, and confirmation at lager 
Nf is required. 



3. Equation of state 

The equation of state (EOS) is one of the most 
basic information in a phenomenological study 
of QGP. In this section, we discuss the latest 
results for EOS from full QCD. Among various 
methods to compute the pressure, the integral 
method Q is commonly used in recent investiga- 
tions. This method is based on the equation, 
p = — /, valid for large homogeneous systems, 
where / — {—T/V)hiZ is the free energy density. 
Because the derivatives of the partition function 
can be expressed by expectation values of oper- 
ators, which are computable by a Monte-Carlo 
simulation, we obtain the pressure by integrat- 
ing this expectation value in the parameter space. 
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Figure 5. Pressure calculated with RG gauge and clover quark actions on a IG'^ x 4 lattice, as a function 
of K (left) and /3 (right). The short lines on the symbols in the right figure express the slope of the 
pressure in the /3-direction, calculated independently ]l7t . 



For the case of the Wilson quark, we have 

(value at T = 0)} (3) 



rp4 rp4 t 



with d^ = {dj3' ^clK') on the integration path. 
The starting point of the integration path should 
be chosen such that p « there. 

3.1. Equation of state with Wilson quarks 

The CP-PACS collaboration carried out a sys- 
tematic calculation of EOS iw Nf = 2 full QCD, 
using the RG gauge and the clover quark actions 
on 16^ X 4 and 16^ x 6 lattices Q . Simulations at 
T = for subtraction and also to fix the scale and 
the line of constant physics, were made on a 16^ 
lattice. The derivatives of the action in (3 and K 
were computed using a U(l) noise method. They 
found that the derivative in the iC-direction gives 
a smaller error in the final pressure. Therefore, 
they chose the integral path in the X-direction. 
Figure ^ (left) is the if-dependence of p/T^ for 
various (3. In order to check the reliability of this 
method, they also show the data as a function of 
/3 in Fig. || (right). The short lines on the symbols 
are the slope of the pressure in terms of /?, which 
are measured directly from {dS/ df3) . We see that 
the slopes well describe the pressure curves which 
are computed from an independent measurement 
of (dS/dK). This means that the value of the 
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Figure 6. mps/Tny constant lines and T/Tpc 
constant lines for an A^^ = 4 lattice in the K) 
parameter space. The values of T/Tpc for the 
dashed lines are given on the right edge of the 
figure0. 



pressure is independent of the choice of the inte- 
gration path in the parameter space. 

The data in this form is not yet useful for phe- 
nomenological applications. We wish to calcu- 
late the temperature dependence for each physi- 
cal system, i.e., on a line of constant physics. The 
CP-PACS collaboration has chosen the condition 
mps/mv — constant, determined at T = 0, to 
define the lines of constant physics. The relation 
between the simulation parameters and the phys- 
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Figure 7. Pressure calculated with the RG gauge Figure 8. Energy density as a function of T/Tpc, 
and clover quark actions on 16'^ x 4 and 16'^ x 6 calculated with RG gauge and clover quark ac- 

lattices, as a function of T/Tpdjl^- tions on 16"^ x 4 and 16"^ x 6 lattices [pT|. 



ical parameters is then given by Fig. ^ The solid 
lines are the lines of constant physics. The dashed 
lines are T/Tpc = constant lines, where Tpc is the 
pseudo-critical temperature on the same line of 
constant physics and the scale was fixed by my- 
The dashed line T /Tpc — 1 itself was determined 
by fitting the pseudo-critical transition points (di- 
amonds) obtained on an 7V( = 4 lattice. 

Interpolating the data in Fig. ^ yields Fig. 
in which the pressure is given as a function of 
temperature, for each line of constant physics. 
Filled (open) symbols are the results of iVt = 4 
(6). Different shapes of the symbol correspond to 
different values of mps/mv, i.e., different quark 
masses. This figure shows that the pressure is 
almost independent of the quark mass in a wide 
range of mps/my- On the other hand, the Nt- 
dependence is sizeable: there exist large lattice 
artifacts in the iVt = 4 results. 

The energy density is obtained by combining 
the results for the pressure and the results for the 
quantity; 



rp4 



1 



(value at T = 0)] . 



dK dS, 



(4) 



Here, the beta functions, and can be 

' ^ da da ^ 

determined by measuring the differentials of (3 
and K in terms of the mya, on each mps/mv — 



constant line. 

Figure || is the final results for the energy den- 
sity as a function of the temperature. They in- 
dicate that the quark mass dependence for the 
energy density is also small. 

Previous results from the staggered quark 
showed a peak just above Tc on an TVj = 4 lat- 
tice! 14 . We do not see such a peak in Fig. || from 
the improved Wilson quark. Also the peak disap- 
pears at A^t = 6 for the staggered fermion, thus, 
the peak seems to be a lattice artifact. 

We find that the magnitude and T-dependence 
of EOS on A^t = 6 lattices are quite similar be- 
tween improved Wilson and staggered quarks. 
The quark mass dependence was also investigated 
with the staggered quark; similar to the case of 
the improved Wilson quark, the difference be- 
tween niq/T = 0.075 and 0.15 is found to be small 
within errors jlj]. 

Short lines on the right axes of Figs. and ^ 
are the Stefan-Boltzmann (SB) values at Nt = 4, 
6 and in the continuum limit. Both A't = 4 and 6 
results do not indicate a clear approach towards 
these values in the high temperature limit. A 
clear deviation from the SB limits at finite Nt 
was also observed in a quenched study of EOS 
with the RG-improved action^. Further study is 
necessary to clarify this phenomenon. However, 
quenched studies show that EOS extrapolated to 
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Figure 9. Pressure from the p4-improved stag- 
gered action on a 16'^ x 4 lattice as a function of 
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Figure 10. Energy density from the p4-improved 
staggered action on a 16"^ x 4 lattice as a function 
of T/Tp,^. 

the continuum limit iVt = c» is consistent with 
the continuum SB value in the high temperature 
limit. From Figs, ^and |[ we note that the results 
for Nf — 6 are already close to the continuum SB 
limit. An accurate continuum extrapolation may 
be possible already with data at A^t ~ 8 or 10. 

3.2. 7V/-dependence of equation of state 

So far, we have discussed simulations with 
2 dynamical quarks, which are identified as 
light up and down quarks. In this subsection, 
We review results on the effect of the strange 
quark. The Bielefeld group compared EOS for 
Nf = 2,3 and 2 + 1 on a 16^ x 4 lattice, us- 
ing the p4-improved staggered quark action with 



the Symanzik improved gauge action|l^. Nt- 
dependence is known to be small for this action 
in the Stefan-Boltzmann limit. They computed 
EOS at rUg/T = 0.4 for Nf = 2 and 3, and at 
mu,d/T = 0.4 and m./T = 1.0 for A^/ = 2 + 1. 

Figure ^ shows the results of pressure. A clear 
A^/-dependence is visible; the pressure becomes 
larger as Nf is increased. 

Clear A^/-dependence was observed also in the 
energy density. They also observed peaks of the 
energy density near Tpc, similar to the result of 
the MILC collaboration at Nt = 4. 

They attempted to estimate the energy den- 
sity in the chiral limit by subtracting the contri- 
bution of {dS/d{mqa)) in the calculation of the 
(e — 3p)/T^, since this term vanishes when the 
quark mass is zero. Quark mass dependence was 
assumed to be small in the pressure. Figure |l^ 
shows the resulting estimate of e/T^ in the chiral 
limit. The peak disappeared in this figure. There- 
fore, in contrast to the case of Fig. ^, noticeable 
quark mass dependence was suggested for the en- 
ergy density. Further study is necessary to clarify 
this point. 



4. Summary 

We reviewed progress of the numerical study of 
finite-temperature QCD with dynamical quarks 
in the past year, focusing on the topics of EOS. 

For the transition temperature, recent results 
for Nf — 2 QCD, computed with improved stag- 
gered and improved Wilson fermion actions on 
Nt = i lattices, agree well with each other. In the 
chiral limit, the transition temperature is about 
170 MeV. This value is higher than the previous 
results from the standard staggered fermion at 
A't = 4 and 6. Confirmation of this result re- 
quires further study at lager Nt- 

The first calculation of EOS with Wilson-type 
fermion was performed for Nf = 2 at Af = 4 and 
6, using an improved Wilson fermion. The results 
of EOS for A^t = 6 are quite similar between the 
improved Wilson and the staggered fermions, and 
are close to the continuum Stephan-Boltzmann 
value at high temperatures. This may suggest 
that a small A'f ~ 8 or 10 is already sufficient 
for a precise calculation of EOS. On the other 
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hand, a clear 7V/-dependence was reported for 
EOS. This imphes the importance of including 
dynamical strange quark in the simulations. 

In summary, we have obtained two keywords 
for the next step of full QCD thermodynamics: 
larger Nt with improved actions towards the con- 
tinuum limit, and simulations with the dynamical 
strange quark. Studies in these directions are now 
starting. 

I would like to thank A. Ukawa and K. Kanaya 
for helpful suggestion and comments on the 
manuscript. I also thank F. Karsch for commu- 
nication on their results and useful discussions. 
This work is supported by JSPS. 
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